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ABSTRACT 

We solve equations, describing in a simplified way the newtonian dynamics of a self- 
gravitating nonrotating spheroidal body after loss of stability. We find that contrac- 
tion to a singularity happens only in a pure spherical collapse, and deviations from the 
spherical symmetry stop the contraction by the stabilising action of nonlinear non- 
spherical oscillations. A real collapse happens after damping of the oscillations due 
to energy losses, shock wave formation or viscosity. Detailed analysis of the nonlinear 
oscillations is performed using a Poincare map construction. Regions of regular and 
chaotic oscillations are localized on this map. 

Key Vifords: gravitation - instabilities. 



1 INTRODUCTION 

Dynamic stability of spherical stars is determined by an av- 
erage adiabatic power 7 — g'i°g"^ \ s- For a density distribu- 
tion p = poip{m/M), the star in the newtonian gravity is sta- 
ble against dynamical collaps e when (7 — |)f .^^^I^ m^ > 
(|Zeldovich fc Novikovlll967l : lBisnovatvi-Koganlll989i ). This 
approximate criterium becomes exact for adiabatic stars 
with constant 7. Here po is a central density, M is a stellar 
mass, m is the mass inside a Lagrangian radius r, so that 
m = 47r pr^dr, M — m{R), i? is a stellar radius. Collapse 
of a spherical star may be stopped only by a stiffening of the 
equation of state, like neutron star formation at late stages 
of evolution, or formation of fully ionized stellar core with 
7 = I at collapse of clouds during star formation. Without 
a stiffening a spherical star in the newtonian theory would 
collapse into a point with an infinite density (singularity). 

In the presence of a rotation a star is becoming more 
dynamically stable against collapse. Due to the more rapid 
increase of a centrifugal force during contraction, in com- 
parison with the newtonian gravitational force, collapse of a 
rotating star will be always stopped at finite density by cen- 
trifugal forces. Here we show, that deviation from the spher- 
ical symmetry in a non-rotating star with zero angular mo- 
mentum leads to a similar stabilization, and non-spherical 
star without dissipative processes never will reach a singu- 
larity. Therefore collapse to a singularity is connected with 
a secular type of instability, even without rotation. 

When a uniform sphere of cold gas is set into free-fall 
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collapse under its self-gravity its shape is unstable. The lin- 
ear instability of large-scale modes (seco nd-order harmonics) 
during the collapse was discovered by iLvnden-Beli (Il964l) 
and f or the ca se of non-rotating sphere bv iLvnden-Belll 
(Il979t ) (see also iLvnden-Belll ()l996h ). Numerical investiga- 
tions of collaps ing pressureless sph eroi ds have been done 
in the papers of|Lxn dcn-Bcii' (| 19641 ) and lLin. Mestel fc "Sh^ 
|l965l ). In iLvnden-Bell (19791) it was also shown that the 
pressure prevents the development of large-scale shape insta- 
bility if i nitially th e gravity is more than three-fifths pressure 
resisted. i Fujimotol (1968 ) derived the equations of motion for 
a rotating ellipsoid from a system of hydrodynamical equa- 
tions. However, the accounting for pressure effects in his 
approach was inconsistent with thermal processes, leading 
to the wrong results for the dynamical behaviour of a sys- 
tem w ith radiative losses (see Bisnovatvi-Kogan fc Tsupkol 
(|2005l )). Using the variational principle and deriving the 
equation for the entropy function from energy balance give 
th e correct relations f or the pressure and the total energy. 
In lRosensteel fc TranI (|l99ll ) the virial equations for rotating 
Riemann ellipsoids of incompressible fluid are demonstrated 
to form a Hamiltonian dynamical system. 

We calculate a dynamical behavior of a non-spherical, 
non-rotating star after its loss of a linear stabil- 
ity, and investigate nonlinear stages of contraction. We 
use approximate system of dynamic equations, describ- 
ing 3 degrees of freedom of a uniform self-gra v itatin g 
compressible ellipsoidal body ([Bisnovatvi-Koganl (|2004l ). 
iBisnovatvi-Kogan fc Tsupkol ( 20051 )). We obtain that the de- 
velopment of instability leads to the formation of a regularly 
or chaotically oscillating body, in which dynamical motion 
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prevents the formation of the singularity. We find regions of SGm °f du 5 1 e 

chaotic and regular pulsations by constructing a Poincare b — — b / , „ — r-- + - — — , , , . ,„ , (9) 

„ , „;..., . . , 2 / (o^ +u)A 3m 6 (aoc)i/'^ 

diagram tor ditterent values ot the mitial eccentricity and p 

initial entropy. For simplicity we restrict ourself by calcu- ^ 

lating only spheroidal figures with 7 = |, and only briefly ^ SGm f du 5 1 e 

represent results for 7 = |. At the end we discuss qualita- ~ 2 I (c^ + it)A ~'~ 3m c {abcy/^ ' 

tively effects of general relativity in a non-spherical collapse 

of a non-rotating body. 2 , 2 , -,n'2 , w 2 , s 

A — (a + u)(b + u)[c + u) . 



2 EQUATIONS OF MOTION 

Let us consider 3-axis ellipsoid with semi-axes a ^ b c: 



, 



(1) 



and uniform density p. A mass m of the uniform ellipsoid is 
written as (1/ is the volume of the ellipsoid) 

47r 

m — pV = — p abc . (2) 

Let us assume a linear dependence of velocities on coordi- 
nates 



3 DIMENSIONLESS EQUATIONS AND 
NUMERICAL RESULTS 

To obtain a numerical solution of equations we write them 
in non-dimensional variables. Let us introduce the variables 

~ t _ jj ^ - 

to' ao ' ao ' CIO ' 

m . P r'r U - Eth - e 

m = , p = U ^ Eth^ -j^, £ = — . 

mo po Uo Uo £0 

The scaling parameters to, ao, mo, po, Uo, £0 are con- 
nected by the following relations 



by 



cz 
c 



(3) 



Th e gravitational energy of the uniform ellipsoid is defined 
as l|Landau fc Lifshitzlll99j ): 



10~ 



du 



y^ia'^ + u)(b'2 + u){c'2 + u) 



(4) 



The equation of state P = Kp'' is considered here, with 
7 = 4/3. Note that the case 7 = 5/3 was considered by 
iBisnovatvi-KoganI l|2004l ) and iBisnovatvi-Kogan fc Tsupkol 
but this case is not interesting for the present 
work, because isentropic spherical star with 7 = 5/3 al- 
ways stops contraction, and never suffers collapse to sin- 
gularity. A spherical star with 7 = 4/3 collapses to singu- 
larity at small enough K, and we show here, how devia- 
tions from a spherical form prevent formation of any singu- 
larity. For 7 = 4/3, the thermal energy of the ellipsoid is 
Eth ~ V"^''^ ~ (afcc)"^''^, and the value 

s = Eth{abc)^^^' = 3[^y''K 

remains constant in time. A Lagrange function of the ellip- 
soid is written as 



L — Ukin — Upot , Upot — Ug + Eth , 



(5) 



to 



Gmo 



, Uo 



ao 



Po 



1 rr 

-, £0 = Uoao- 



System of non-dimensional equations: 



du 5 1 £ 

+ 



3m 

2~"7 (a2 + M)A ' 3m a (afec) 1/3 




" 3m f du 5 1 £ 

" 2" J (62 + u)A ^ 3^ 6 (a6c)i/3 ' 



du 5 1 £ 

+ 



3m 

2^""] (c2 -I- ' Ikn c [abcy/^ ' 





(11) 



(12) 



(13) 



(14) 



= [a^ + u){b^ + u){^ + u) . 

In equations (I12|l - (|14| ) only non-dimensional variables 
are used, and "tilde" sign is omitted for simplicity in this sec- 
tion. The non-dimensional Hamiltonian (or non-dimensional 
total energy) is: 

H = [/fc,„ + Ug+Eth^'^ {a +1?+ c")- 
10 



3m2 

To" 



du 



+ 



^(a2-f u)(62-f m)(c2-)-u) {abcy/'i 



(15) 



Ukin ^\p i {vl + vl + vl)dV ^ ^ (d^ + b^ + c) , (6) 



Eth = 



(afoc) 1/3 



(7) 



Equations of motion describing behavior of 3 semiaxes 
(a, b, c) is obtained from the Lagrange function (O in the 
form 



3Gm / du 
a = — a 



+ 



5 1 



2 J (a2 + u)A 3m a (abcy/^ 




(8) 



In case of the sphere {a — b — c, d — b — c) the non- 
dimensional Hamiltonian and non-dimensional equations of 
motion reduce to: 



Tr 3 .2 3/2 5 

^=10'"" -5^r "3"'' 



a = 3_ fm- _ £ 

ma2 V 3 



(16) 



(17) 



As follows from p6|l . (|17p for the given mass there is only 
one equilibrium value of e 



3m^ 



(18) 
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at which the spherical star has zero total energy, and it may 
have an arbitrary radius. For smaller e < Eeq the sphere 
should contract to singularity, and for e > Eeq there will be 
a total disruption of the star with an expansion to infinity. 
We solve here numerically the equations of motion for a 
spheroid with a = b ^ c, which, using l|12[) - l|14| l. are written 
for the oblate spheroid with k = c/a < 1 as 



.. 3 

a = - 



2a2(l-fc2) 
a2(l--fc2) 



arccos k 
k arccos k 



+ 



5 1 £ 

3m a [a'^cY^^ ' 
5 1 e 



+ 



3m c (a'^c)^/^ ' 



and for the prolate spheroid k = c/a> 1 as 



2a2(fc2 



k- 



cosh ^ k 



Vfc2 - 1 
k cosh~^ k 



+ 



+ 



5 1 £ 

3m a (a2c)i/3^ 
5 1 £ 



'a2(fc2-l)[" Vfc2 - 1 J ' 3m c (a2c)i/3- 

It is convenient to introduce variables 
5 £ 



£* 



3 m2 



tJm. 



(19) 
(20) 

(21) 
(22) 

(23) 



In these variables the equations (I19|l - p2|l . p7|) are written 
as (omitting subscript "*") 



2a2(l-fc2) 
3 



a2(l 



k 



1 - 



arccos k 
VI - fc2 

k arccos k 
Vl - fc2 



+ 



+ 



a (a2c)V3' 
1 £ 



for the oblate spheroid k = c/a < 1, 



2a2(fc2 
3 



cosh ^ fc 



a2(fc2 



1 



fc cosh"^ fc 



+ 



c (a2c)i/3 



1 



(24) 
(25) 

(26) 
(27) 

(28) 

for the sphere, where the equilibrium corresponds to e^q — 1. 
Near the spherical shape we should use expansions around 
k — 1, what leads to equations of motion valid for both 
oblate and prolate cases 

k 



VP" 



+ 



a (a2c)i/3 

1 £ 
c (a2c)V3 



for the prolate spheroid k = c/a > 1, and 

1"£ 



1 



£ 3 
3 + 5 



1 



/4£ _ 4\ ]_ 
V 3 ~ 5J . 



+ 



In these variables the total energy is written as 
m^ 

and omitting we have 
3 arccos k 



(29) 



(30) 
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3 

+ 



H 



a 

"5" 



c 

10 



5a V]~P ■ 5 (a2c)i/3' 

£ 



3 cosh ^ fc ^ 3 



5a VF 



1 



5 (a2c) 



1/3 ' 



(oblate) 
(prolate) 




Figure 1. Example of regular motion of spheroid with 7 = 4/3, 
H = —1/5, e = 2/3. This motion corresponds to full line on the 
Poincare map in Fig. 4. 



i?=^-^(l-£). 



10 

•2 



5a 
• 2 



H 



a c 
"5" + 10 



_3_ 

5a 



(sphere) 



5_ 2S^ 
3+15 



(31) 



+ ■ 



3£ 

5a 



S 2S^ 
3 + - 



1 — k, (around the sphere) 



< 1. 



Solution of the system of equations (|19p - (|22p was performed 
for initial conditions at t = 0: co = 0, diff'erent values of 
initial ag, do, fco, and different values of the constant pa- 
rameter £. Evidently, at fco = 1, do = 0, £ < 1 we have the 
spherical collapse to singularity. The most interesting result 
was obtained at fco 7^ 1, and all other cases with deviations 
from spherical symmetry. No singularity was reached in this 
case for any £ > 0. It is clear, that at £ = a weak singu- 
larity is reached during formation of a pancake with infinite 
volume density, and finite gravitational force. At £ > the 
behavior depends on the value of the total energy H: at 
H > we obtain a full disruption of the body, and at H < 
the oscillatory regime is established at any value of £ < 1. At 
£ ^ 1 the total energy of spheroid is positive, H > 0, deter- 
mining the full disruption. The case with H — is described 
separately below. 

At H < the type of oscillatory regime depends on 
initial conditions, and may be represented either by regu- 
lar periodic oscillations, or by chaotic behavior. Examples 
of two types of such oscillations are represented in Figs. 1-2 
(regular, periodic), and in Fig. 3 (chaotic). For rigorous sep- 
aration between these k inds of oscillations we use a met hod 
developed by Poincare l|Lichtenberg fc Liebermanlll983l ). 

3.1 The case H=0 

In the case of a sphere with zero initial radial velocity we 
have equilibrium state with e = Eeq ~ 1. Let us consider 



the case of d 7^ and £ < 1. It follows from (|3H) . that for 
non-zero initial velocity d entropy function e is strictly less 
than unity. The fate of such sphere depends on the sigh of d. 
At d < the star collapses to singularity, and at d > there 
total disruption happens, with zero velocity at infinity. 
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Figure 2. Example of regular motion of spheroid with 7 = 4/3, 
H = —1/5, e = 2/3. This motion corresponds to the point inside 
the regular region on the Poincare map in Fig. 4. 




Figure 3. Example of chaotic motion of spheroid with 7 = 4/3, 
H = —1/5, e = 2/3. This motion corresponds to gray points on 
the Poincare map in Fig. 4. 



For the initial spheroid with a = b c, and zero ini- 
tial velocities d = c = the entropy function is also less 
than 1. In this case the value of the entropy function e is 
uniquely determined by the deformation 5 = {a — c)/a. This 
dependence can be found explicitly from (|31|l : for small S we 
obtain 



e = 1- 



45 



(32) 



Thus even in the case of zero initial velocities, e is less than 
unity for a spheroid, and reaches unity only for the sphere. 
For the same e the deformation 5 of the zero energy body at 
rest have two states with S — ±§-\/5(l — e), corresponding 
to oblate and prolate spheroids. If we set non-zero initial 
velocities in case of spheroid, entropy function will be even 
less. Calculation of the motion in this case leads finally to the 
expansion of the body, with oscillating behavior of S around 



zero (oblate - prolate oscillations). It takes place both under 
the condition of the initial contraction or initial expansion. 



4 THE POINCARE SECTION 

To investigate regular and c haotic dynamics we use 
the m ethod of Poincare section l|Lichtenberg fc Liebermanl 
I1983D and obtain the Poincare map for different values of 
the total energy H. Let us consider a spheroid with semi- 
axes a = b ^ c. This system has two degrees of freedom. 
Therefore in this case the phase space is four-dimensional: 
a, a, c, c. If we choose a value of the Hamiltonian Ho, we fix 
a three-dimensional energy surface H{a,a, c, c) — Ho- Dur- 
ing the integration of the equations (|19p - (|22p which preserve 
the constant H, we fix moments ti, when c = 0. At these 
moments there are only two independent values (i.g. a and 
d), because the value of c is determined uniquely from the 
relation for the hamiltonian at constant H. At each moment 
ti we put a dot on the plane (a, a). 

For the same values of H and e we solve equations of 
motion (|19p - (|22p at initial c = 0, and different a, a. For 
each integration we put the points on the plane (a, a) at the 
moments ti. These points are the intersection points of the 
trajectories on the three-dimensional energy surface with a 
two-dimensional plane c = 0, called the Poincare section. 

For each fixed combination of e, we get the Poincare 
map, represented in Figs. 3-6. Condition c = splits in two 
cases, of a minimum and of a maximum of c. The Poincare 
maps are drawing separately, either for the minimum, or for 
the maximum of c, and both maps lead to identical results. 
The regular oscillations are represented by closed lines on 
the Poincare map, and chaotic behavior fills regions of fi- 
nite square with dots. These regions are separated from the 
regions of the regular oscillations by separatrix line. 



4.1 The bounding curve 

Actually, the variables a and a cannot occupy the whole 
plane (a,d) : < a < 00, —00 < d < +00. Let us obtain a 
curve bounding the area of the values a and d. Let a function 
$(a,d, c) in the variables (^5)1 . (I5U|| be 



$(a,d,c) = — {(i^ + b^ + c^) 



du 



2,e 



10 y J (a^ + u)aP + u)(c^ + u) 5(abc)i/3 




H (33) 



aX a — b, a = b, c = 0, and fixed value of e and H. The 
equation for the bounding curve /(a, d) = at given e and 
H is determined from the following system of equation 



d 

<E>(a, d, c) = , ^<I>(a,d, c)=0. 

In case of c < a this system of equations has a form 

3e 



1.2 3 arccos(c/a) 
5" "5 ^ 5(a2c)i/3 

3 arccos(c/a) 3 
-c— ^hri^ -f - ■ 



H = Q, 



1 1 



3e 



5 (a2 - c2)3/2 ' 5 a2 



3 c5(a2c)i/3 



= 0. 



(34) 

(35) 
(36) 
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Figure 4. The Poincare map for five regular and two chaotic tra- 
jectories in case of 7 = 4/3, H = —1/5, e = 2/3. The (a, a) values 
are taken in the minimum of c. Full black line is the bounding 
curve. The point inside the regular region corresponds to coherent 
oscillations with the same period for a and c values, represented 
in Fig.2. 



Figure 6. The Poincare map for five regular and two chaotic tra- 
jectories in case of 7 = 4/3, H = —1/5, £ = 2/3. The (a, a) values 
are taken in the maximum of c. Full black line is the bounding 
curve. The point inside the regular region corresponds to coherent 
oscillations with the same period for a and c values, represented 
in Fig.2. 
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Figure 5. Zoom of previous figure. 



This system is solved numerically. The second equation 
doesn't depend on a. We set a and obtain corresponding 
value of c from the second equation. Then we substitute a 
and c into the first equation and find a. Thus we obtain the 
point {a, a). Changing a, we obtain numerically the curve 
/(a, d) = 0. This bounding curve is painted in Figs. 4-8 by a 
heavy line. 



5 DISCUSSION 

The main result following from our calculations is the indi- 
cation to a degenerate nature of formation of a singularity 
in unstable newtonian self-gravitating gaseous bodies. Only 
pure spherical models can collapse to singularity, but any 
kind of nonsphericity leads to nonlinear stabilization of the 
collapse by a dynamic motion, and formation of regularly 
or chaotically oscillating body. This conclusion is valid for 
all unstable equations of state, namely, for adiabatic with 



Figure 7. The Poincare map for two chaotic trajectories in case 
of 7 = 4/3, H = —1/2, e = 1/6. The (a, d) values are taken in 
the minimum of c. Full black line is the bounding curve. 



7 < 4/3. In addition to the case with 7 — 4/3, we have 
calculated the dynamics of the model with 7 = 6/5, and 
have obtained similar results. The Poincare map for this 
case is represented in Fig. 9. For 7 = 6/5 we have the entropy 
function e = Eth{abcf'^^ , nondimensional entropy function 
e = e/Uoa^^^ , see (|lip . For the spherical star with 7 = 6/5 
the nondimensional Hamiltonian instead of (1161) is 



H 



3m , 

c 

10 



So £ 
^■m + -375. 

5a a-^'-^ 



(37) 



Note, that region of chaotic behavior on the Poincare map 
is gradually increasing for 7 = 4/3 with decreasing of 
the entropy e and the total energy H. At e = 9/10 
and H — —3/50 we have found only regular oscillations 
(Fig.8), at £ = 2/3 and H = -1/5 both kind of oscil- 
lations are present (Figs. 4-6), and only chaotic behavior 
is found at e — 1/6 and H = —1/2 (Fig.7). We connect 
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Figure 8. The Poincare map for six regular trajectories in the 
case of 7 = 4/3, H = —3/50, e = 9/10. The (a, d) values are 
taken in the minimum of c. Full black line is the bounding curve. 
The point inside the regular region corresponds to coherent oscil- 
lations with the same period for a and c values, similar to those 
represented in Fig. 2. 



ijBisnovatvi-Kogan fc Tsupkd l2005h this chaotic behavior 
with development of anisotropic instability, when radial ve- 
locities str ongly exceed the transversal ones (|Antonovlll973 : 
lFridman~fc Polyachonko 1985). 

In reality a presence of dissipation leads to damping of 
these oscillations, and to final collapse of nonrotating model, 
when total energy of the body is negative. In the case of core- 
collapse supernova the main dissipation is due to emission of 
neutrino. The time of the neutrino losses is much larger than 
the characteristic time of the collapse, so we may expect that 
the collapse leads to formation of a neutron star where non- 
spherical modes are excited and exist during several seconds 
after the collapse. In addition to the damping due to neu- 
trino emission the shock waves will be generated, determin- 
ing highly variable energy losses during the oscillations. Be- 
sides viscosity and radiation which can damp the ellipsoid- 
like motions and allow collapse, there is the possibility that 
inertial interactions with higher order modes that must be 
present in the real stratified bodies may cause an inertial 
cascade and drain the energy from the second order modes 
an a non-secular way that does not depend on dissipative co- 
efficients. It is very seductive to connect chaotic oscillations 
with highly variable emission observed in the p rompt gamma 
ray e mission of cosmic gamma ray bursts ([Mazets et al.l 
Il98ll ). A presence of rotation and magnetic field strongly 
comp licate the picture of the core-coll apse supernova explo- 
sion l|Ardelian. Bisnovatvi-Kogan fc M ois eenko. 2005 ). 

In this paper we consider in details only spheroidal 
bodies. In reality the spheroid will become a triaxial el- 
lipsoid during the motion. In addition to spheroids we 
calculate many variants w i th tr iaxial figures (see also 
iBisnovatvi-Kogan fc Tsupkd (|2005l )). Qualitatively we ob- 
tain the same results for ellipsoids: no singularity was 
reached for any e > and establishing oscillatory (regular 
or chaotic) regime under negative total energy prevents the 
collapse. However, in the case of the ellipsoid with semi-axes 




Figure 9. The Poincare map for one chaotic and four regular 
trajectories in case of 7 = 6/5, m = 1, H = —3/50, e = 27/50, 
see I I37I I. The (a, d) values are taken in the minimum of c. Full 
black line is the bounding curve. The point inside the regular 
region corresponds to coherent oscillations with the same period 
for a and c values, similar to those represented in Fig. 2. 

a 7^ 6 7^ c we have a system with three degrees of freedom 
and six-dimensional phase space. Therefore we could not 
carry out rigorous investigation of regular and chaotic types 
of motion by the constructing Poincare map as it was done 
for spheroid with two degrees of freedom, and restricted our- 
selfs by a description of the spheroidal case. 

In the frame of a general relativity dynamic stabilization 
against collapse by nonlinear nonspherical oscillations can- 
not be universal. When the size of the body approaches grav- 
itational radius no stabilization is possible at any 7. Nev- 
ertheless, the nonlinear stabilization may happen at larger 
radii, so after damping of the oscillations the star would col- 
lapse to the black hole. Due to development of nonspherical 
oscillations there is a possibility for emission of gravitational 
waves during the collapse of nonrotating stars with the in- 
tensity similar to rotating bodies, or even larger. 

Account of general relativity will introduce a new non- 
dimensional parameter, which can be written as Pq — — 2 — — . 
The fate of the gravitating body will depend on the value 
of this parameter, and we may expect a direct relativistic 
collapse to a black hole at increasing pg , approaching unity. 
It is known, that a nonrotating b lack hole is characterized 
only by its mass (|Zeldovich fc Novikov 1967). In absence of 
other dissipative processes, the excess of energy, connected 
with a nonspherical motion should be emitted by gravita- 
tional waves during a formation of the black hole. 
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